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After many years of speculation, recent observations have confirmed the association of gamma- 
ray bursts with core-collapse supernova explosions from massive stars. This association carries 
with it important consequences. The burst relativistic jet has to propagate through the cold dense 
stellar material before it reaches the transparency radius and the burst photons are produced. This 
propagation is likely to affect the initial properties of the jet, shaping it and changing its energy 
composition. The variability injected at the base of the jet is also likely to be erased by the jet-star 
interaction. Despite this, GRBs seem to have remarkably predictable properties once the radiative 
phase sets in, as emphasized by the recent discovery of several tight correlation between spectral, 
geometric and energetic properties of the jet. In this contribution we discuss the jet interaction 
with the star, emphasizing its time-dependent properties and the resulting energy distribution. We 
finally emphasize the surprising predictability of jet and radiation properties outside the star and 
underline its implication for standardizing the GRB candle. 


I. INTRODUCTION 

After many years of speculation Hi , the associa¬ 
tion of (at least some) long gamma-ray bursts (GRBs) 
with supernova (SN) explosions has been established. 
The first association was between GRB 980425 and 
SN1998bw [ 3 . It was a problematic and highly de¬ 
bated association since it implied a highly peculiar na¬ 
ture for GRB 980425: an event with a normal fluence 
but a distance 100 times smaller then other GRBs, 
and therefore some four orders of magnitude fainter. 
In addition, GRB 980425 did not have a normal af¬ 
terglow in any band. Nevertheless, the peculiar na¬ 
ture of the SN explosion and the temporal and spatial 
coincidence argued strongly for the association. Af¬ 
ter GRB 980425/SN1998bw, circumstantial evidence 
for more associations was claimed for higher redshift 
events, based on the appearance of a red bump in the 
late afterglow light curve 0. 

More recently two additional events have strength¬ 
ened the association, including - at least in one case 
- one GRB that has all the characteristics of classical 
cosmological events. GRB 030329 was the brightest 
GRB detected by HETE-2 Q. It had a small redshift 
2 : = 0.168 and multi-epoch afterglow spectroscopy 
revealed the emergence of a SN spectrum ^proxi- 
mately 15 days after the burst explosion The 

spectrum and the overall supernova properties were 
remarkably similar to those of SN1998bw. Finally, 
GRB 031203, detected by INTEGRAL, looks like a 
twin to GRB 980425/SN1998bw 

The association bears new questions and riddles 
with it. First we know that the 7 -ray photons are 
produced by ultra-relativistic material. The progeni¬ 
tor star is however dense and massive, and any outflow 
collecting matter geometrically from the star will be 
slowed-down to sub-relativistic speed very effectively. 


The entrained mass can be at most a small fraction 
of a solar mass Mq < 5.5 x 10 ~^ E^i F^^ Mq. Simula¬ 
tions show that such a small contamination is possible 
(see below). The jet that reaches the surface of the 
star is then shaped by the interaction with the star 
itself. One would expect that, depending on the type 
of stellar progenitor, jet initial conditions and ener¬ 
getics, a great diversity of jets would emerge. After 
all type Ibc and type II supernovae show a remarkable 
diversity. 

Even though at first sight GRB observations are 
characterized by a huge diversity, remarkable correla¬ 
tions have been discovered among jet properties. The 
first such correlation was discovered between the jet 
7 -ray energy output (not corrected for the beam ing ) 
and the break time in the afterglow light curve [inLlll| . 
This correlation can be understood in terms of jets 
with the same total energy but beamed into different 
opening angles terms of structured jets with 

an energy distribution such that dE/dVL oc 9 -2 [12. 

A further correlation was discovered between the 
typical frequency of photons in the prompt emission 
(the peak of the vF{v) spectrum) and the isotropic 
equivalent energy output in the prompt phase p^ . 
According to this correlation the brighter the burst, 
the larger the typical frequency of photons. Einally, 
the tightest correlation of all was recently discovered 
between the typical photon energy and the beaming 
corrected 7 -ray output H (under the assumption of 
a uniform jet). All these correlations, and especially 
those involving the photon frequency, have not been 
understood so far. 

In this contribution we discuss the jet propagation 
inside the progenitor star |15|. We consider in par¬ 
ticular the effect of the diminishing collimating power 
of the star, which is blown apart by the energy dis¬ 
sipated by the jet as it opens its way out of the pro- 
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genitor. We show that this creates time dependent 
effects that may be responsible for the diversity in 
the 7 -ray phase of the burst. Such effects, however, 
occur on short timescales, and therefore are not rele¬ 
vant in the afterglow phase, which is therefore more 
standardized. Finally, we discuss the opportunity of¬ 
fered by the Ghirlanda correlation H to use GRBs as 
standard candles to constrain the expansion history of 
the universe at redshifts much larger than the range 
sampled by type la SNe [T^ . 


II. THE STAR-COCOON-JET 
INTERACTIONS 

Consider a jet generated at the core of a massive 
star (at a radius Ref), characterized by a luminos¬ 
ity Lj, an initial Lorentz factor Fq ff 1, a dimen¬ 
sionless entropy rj = Lj/[M(?) (where M is the 
rest mass ejection rate) and an initial opening an¬ 
gle 00 (or, alternatively, a stagnation pressure Pstag = 
L^Tq/{A-k cRqOq)). The jet attempts to propagate 
through the stellar material. As soon as the head 
of the jet reaches supersonic velocity, a double shock 
structure develops at the head of the jet. The speed 
of the head of the jet (the contact discontinuity) can 
be computed by balancing the pressure of the jet ma¬ 
terial with that of the cold stellar matter. Pressure 
balance reads miH: 


Pi (r/3)f/i -k Pj = p* K{TI3)1 c^ + p* (1) 

where p is the density, h the enthalpy and p the pres¬ 
sure. Subscripts are j for jet material, jt for the jet 
material relative to the head of the jet, h for the jet 
head relative to the star and * for the stellar mate¬ 
rial. This equation, if we neglect the pressure terms 
on both sides, can be solved to 
jet head in the stellar material 


give me morion or rne 

ll3: 


/3/t = 



( 2 ) 


where C = Pj hj/p,,. 

The actual structure of the head of the jet is made 
by a contact discontinuity that separates jet material 
from stellar material, two thin regions of shocked ma¬ 
terial (shocked jet and shocked stellar material) and 
two shocks (a forward bow shock propagating into the 
star and a reverse shock into the jet). The shocked 
material is over-pressured and flows to the sides of 
the jet, feeding a cocoon. 

The cocoon material is over-pressured with respect 
to both the star and the un-shocked jet. It therefore 
expands, driving a shock into the cold stellar mate¬ 
rial and collimating the jet until pressure balance is 
reached. We assume an adiabatic flow, since the jet is 


optically thick inside the star. As a consequence, the 
Lorentz factor scales as F^ oc , where is the 
area of the jet. This relation, which is valid only if 
dissipation is negligible, allows us to compute the jet 
pressure as a function only of Sj and to derive conse¬ 
quently the Lorentz factor. The cocoon evolution can 
be computed by applying the first law of thermody¬ 
namics: 


d{ec Vc) = dQ -pc dVc (3) 

where Cc is the cocoon energy density and Vc is the 
volume of the cocoon. We assume that the cocoon 
pressure is uniform, so that Pc only depends on time. 
The term dQ = Lj{l — fdh)dt is the energy injected in 
the form of shocked jet material. The evolution of the 
cocoon volume 14 reads: 

dV 

-—L = 2Tr R±{r,t)vshir,t)dr (4) 

™ JRo 

where R± is the transverse size of the cocoon (which 
is a function of radius and time) and Vsh is the 
speed of the shock driven by the cocoon pressure into 
the stellar material. It can be computed under the 
Kompaneets approximation [T^ by balancing the co¬ 
coon pressure against the ram pressure exerted by 
the cold stellar material on the expanding cocoon: 
Vsh = ^Cc/ip*- 

The above set of equations form a solvable system 
of differential equations. An analytic solution is not 
possible, especially if the time-dependent evolution of 
the cocoon needs to be taken into account. 

To show some examples of the cocoon evolution we 
assume a star with mass M* = ISM©, with radius 
ii* = 10^^ cm and p* oc r““*. Figuren]shows the co¬ 
coon evolution at fixed times (indicated in the panels). 
The star has a power-law density profile p* oc r~^'^ be¬ 
tween 10® and 10^^ cm. No funnel pre-evacuation is 
assumed. The jet is injected with Lj = 10®^ erg s“^ 
at ro = 10® cm with 9q = 20° and Fq = 1.5. A recolli- 
mation shock immediately collimates and slows down 
the jet (upper left panel of Fig. P). The jet emerges 
on the stellar surface after ^ 5 s. A total energy of 
1.5x10®^ erg has been dissipated to create the jet path 
into the star. Of this energy, ~ 4 x 10®® erg are stored 
into the cocoon plasma, while ~ 10®^ erg are given 
to the star through the cocoon shock. Such energy 
is of the order of the explosion energy of Ic SNe and 
enough to unbind the star. The jet reaches the star 
surface with an opening angle of ~ 1° and a Lorentz 
factor Fbr — 24 and is therefore in causal contact. 

Figure shows a similar simulation but where an 
exponential cutoff of the stellar density of the form 
exp[—r/(10^®cm)] has been added to simulate a more 
realistic profile. The main effect of this change is on 
the shape of the cocoon that opens as the exponential 
cutoff is reached. Also, the jet is somewhat larger 
and has a larger Lorentz factor so that it turns out 
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to be only in marginal causal contact. These results 
compare positively with analytic estimates M 


III. PRECURSORS 

As the jet pierces the stellar surface, it simultane¬ 
ously opens a channel for the release of the cocoon. 
The cocoon material therefore expands out of the star. 
The actual dynamics and radiative properties of the 
cocoon are somewhat uncertain and several different 
predictions have been made 113, m mm. what 
is clear is that a total energy of severalxlO®° ergs is 
released almost isotropically. 

Observationally, associating a precursor with a 
transient event is challe nging, if at all possible. At¬ 
tempts have been made |24l 125^ adopting conservative 
definitions, and therefore only the tip of the iceberg 
was likely discovered. Figure |3| shows a sub-sample of 
the precursors discovered by Lazzati in an anal¬ 
ysis of ^ 100 bright long BATSE GRBs. He defined 
precursors as activity that is detected before the GRB 
trigger from the same direction as the GRB. The de¬ 
tected precursors are somewhat puzzling. A precursor 
related to the jet release or to any event that takes 
place before the 7 -rays of the prompt phase are re¬ 
leased is supposed to precede the GRB by: 

AT = ^ ~ 0.03 s (5) 

to be compared to the several tens of seconds of de¬ 
lays found observationally. The only wav|^ to under¬ 
stand these precursors as related to the cocoon break¬ 
out is to assume that the jet is initially dark in 7 -rays. 
We argue in the following that the jet spreads during 
the ^ 100 s of the energy release. If the observer is not 
lying on the jet axis, there will be a time interval dur¬ 
ing which there will be no 7 -ray emission, even though 
the precursor had been seen. Once the jet has spread, 
it enters the line of sight to the observer and becomes 
visible, with a delay that depends on the off-axis angle 
of the observer. 


IV. JET BREAKOUT AND SUBSEQUENT 
EVOLUTION 


expands, Cc the cocoon energy density and Cs = c/^/S 
is the sound speed of the relativistic gas. Writing the 
cocoon volume as 14 ~ EcT*, we can obtain the pres¬ 
sure evolution: 


Pc = Pc,hr exp 



( 6 ) 


where Pc,hr is the cocoon pressure at the moment of 
shock breakout. As the cocoon pressure decreases, 
fresh jet material passing through the cocoon is less 
tightly collimated. Under isentropic conditions the 
jet Lorentz factor increases linearly with the opening 
angle and pressure balance yields Oj oc pc imply¬ 
ing an exponentially increasing opening angle of the 

formim 


0j = 0j,hr exp 



(7) 


Dissipative jet propagation gives similar results (with 
merely a different numerical coefficient ~ 0 ( 1 ) inside 
the exponential), provided that T varies roughly as a 
power of Ej. 


V. JET STRUCTURE AT LARGE RADII 


The angular distribution of integrated energy, as 
observed in the afterglow phase, is computed by in¬ 
tegrating the instantaneous luminosity per unit solid 
angle from the moment the jet becomes visible along 
a given line of sight (b.o.s.) until the end of the burst: 


dE 

dD 


/■Tgrb ^2, 

/ dt ~ 

' ^l.o.s. 


dLl 



m 

7r6»2(t) 


dt. 


( 8 ) 


ti.o.s. is obtained by inverting eq.|7| Such integration is 
valid provided that the jet opening angle at time Tqrb 
is smaller than the natural opening angle of the jet: 
Tgrb < 4V3(r*/c) log(6>o/0br). For the fiducial num- 
bers assumed (r* = 10^^ cm, 0ybr = 1° and Oq = 30°) 
this corresponds to ~ 100 s comoving burst duration. 
Assuming a jet with constant luminosity L and for 
all the lines of sight that satisfy tbr < b.o.s. "C Tqrb, 
eq. E gives the jet structure 


As the jet reaches the stellar surface, it clears a 
channel for the cocoon. The cocoon material is there¬ 
fore now free to expand out of the star and its pres¬ 
sure drops. We assume that from this moment on the 
shock between the cocoon and the cold stellar mate¬ 
rial stalls, as a consequence of the dropping cocoon 
pressure. This is equivalent to assuming a constant 
volume of the cocoon cavity inside the star. The pres¬ 
sure drop for the relativistic cocoon can be derived 
through dEc = —CcEcCa dt where Ec is the area of 
the free surface through which the cocoon material 


dE 

dLl 


2V3L 


r* 


9- 


9 


TTC 


7 ,br <0<9o 


(9) 


and dE/dfl ^constant inside the core radius 0j,br- 
This angular dependence, which characterizes a 
“structured jet” or “universal jet” d m m m, 
is of high theoretical interest. Jets with this beam 
pattern reproduce afterglow observations. If the jet 
is powered by fall-back of material from the star to 
the accretion disk, the mass accretion rate would be 
anti-correlated with the radius of the star, for a given 
stellar mass. 
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FIG. 1: Stills from a movie showing the jet and cocoon evolution inside a 15 Mq star with a density profile p* oc r~ 
and a radius of 10^^ cm. The jet is injected at a radius ro = 10® cm with an opening angle 9o = 20°, a Lorentz factor 
To = 1.5 and a luminosity Lj — 10®^ erg s“^. Note that the scale of figures is enlarged to emphasize the shape of the 
cocoon. 
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FIG. 2: Same as Fig. 0 but for a star with an exponential cutoff in the density at large radii. 
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FIG. 3: Atlas of a fraction of the precursors detected in a sample of bright long BATSE GRBs [^ . 


In the above equations we have assumed for simplic¬ 
ity that the jet reaching the surface of the star is uni¬ 
form. As shown by simulations |^. it is more likely 
that a Gaussian jet emerges from the star. On the 
contrary, boundary layers may be produced by the in¬ 
teraction of the jet with the collimating star, resulting 
in edge brightened jets (see the inset of Fig. El). It can 


be shown easily that the 9~^ pattern does not depend 
on the assumption of uniformity. Fig.Elshows the inte¬ 
grated energy distribution for uniform, Gaussian and 
hollow intrinsic jets. Even though small differences 
are present at the edges (the jet core and the out¬ 
skirts), the general behavior is always dE/dLt oc 9~^. 
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FIG. 4: Energy distribution for the afterglow phase for 
three instantaneous beam patterns (see inset). In all three 
cases a well defined dE/dfl oc 9~^ section is clearly visi¬ 
ble. Only the edge of the jet and its core show marginal 
differences. The results shown are for a jet/star with 
L = 10®^ erg s“^, Tqrb = 40 s and r* = 10^^ cm. Inset: 
Instantaneous beam patterns that reach the surface of the 
star. The solid line shows a uniform jet, dashed line shows 
a Gaussian energy distribution, while the dash-dotted line 
shows an edge brightened (or hollow) jet. 


VI. THE E^-Cpeak CORRELATION 

Correlations between observed (and intrinsic) quan¬ 
tities are one of the most remarkable results on the 
general properties of GRB jets. They tell us that, 
despite the varied behavior of the observed jet prop¬ 
erties, they are not entirely random but follow a well 
defined path. Among such correlations 
the more recent and tightest is the so-called Ghirlanda 
correlation 0 that finds a correspondence between 
the beaming corrected y-ray energy release of the 
prompt emission and the typical photon frequency 
(the peak of the vF{y) spectrum). This correlation 
is so tight that, within the present accuracy, it is con¬ 
sistent with being an exact relation. The debate is 
however alive on this issue M- 

The correlation is shown in Fig. [3 On the right of 
the same graph the Amati correlation is shown 
for comparison. The origin of these correlation is not 
understood so far. It is believed that they may be due 
to a less dominant role of synchrotron in the prompt 
radiative phase of the burst [^. 


VII. STANDARD CANDLES AND 
COSMOGRAPHY 

The Ghirlanda correlation is so tight that it can be 
used to standardize the GRB candle M If one is 
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FIG. 5: Rest frame peak energy Ep,,ak = + z) ver¬ 

sus bolometric energy for the sample of GRBs with mea¬ 
sured redshift. Filled circles: isotropic energy corrected 
for the collimation angle by the factor (1 — cos9), for the 
events for which a jet break in the light curve was ob¬ 
served. Grey symbols corresponds to lower/upper limits. 
The Solid line represents the best fit to the correlation, i.e. 
Epeak ~ 480 (E^./lO^^erg)'^’^ keV. Open circles: isotropic 
equivalent energy E~,^iso- The Dashed line is the best £t 
to these points. 


able to measure the break time of the GRB afterglow 
(and therefore constrain the jet opening angle) and the 
typical frequency of photons in the y-ray spectrum, it 
is possible to predict the intrinsic energy radiated in 
y-rays in the prompt phase. This allows one to mea¬ 
sure the luminosity distance independently of redshift. 
If, in addition, a redshift measurement is available, we 
can draw a Hubble diagram and fit cosmological mod¬ 
els to it. 

FigureElshows the GRBs and SNe Hubble diagram 
before and after the empirical corrections (stretch- 
luminosity for SNe and Ghirlanda for GRBs). Fig¬ 
ure [71 shows instead the constraints that can be ob¬ 
tained in the LIm — Ha plane with SNe, GRBs and a 
simultaneous fit of the two samples . Even though 
the simultaneous fit is clearly dominated by SNe, the 
perspective for GRBs, should the Ghirlanda correla¬ 
tion be confirmed, are very good. Their appeal is, 
more than to give an independent confirmation of the 
SNe result, to allow for an expansion to high redshift 
of the universe that can be sampled. This is an ideal 
set-up to investigate on the possible evolution of the 
properties of Ha to understand better the origin of 
dark energy. 

With the launch of Swift, it is envisaged that the 
sample of GRBs with all the necessary data will soon 
become comparable to the SN sample. Fig. 0 shows 
a prediction of how the cosmological constraints from 
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FIG. 6: Classical Hubble diagram in the form of 
luminosity-distance Dl vs redshift ^ for Supernova la and 
GRB. In the top panel the SN la and GRBs are treated 
as standard candles (no corrections applied); for GRBs 
E-y = 10®^ erg is assumed. In the bottom panel we have 
applied the stretching-luminosity and the E^-Ep,,ak re¬ 
lations to SN la and GRBs, respectively. Note that, for 
GRBs, the applied correction depends upon the specihc as¬ 
sumed cosmology, here for simplicity we assume the stan¬ 
dard Q.M = 0.3, Ha = 0.7 cosmology. Both panels also 
show different Dl{z) curves, as labeled. 


GRBs in the Swift era may look. 


VIII. DISCUSSION 

We have discussed the properties of GRB jets 
within the framework of the GRB-SN association. We 
showed that the interaction of the jet with the star 
produces a cocoon that modifies the jet properties. 
Even though the jet may be released in a stationary 
process, with constant properties in time, the inter¬ 
action with the star and the cocoon creates a time- 
dependent opening angle of the jet. In the prompt 
phase the main consequence is that the beginning of 
the 7 -ray emission depends on the observer. In the 



FIG. 7: Gonstraints in the Hm — Ha plane derived from our 
GRB sample (15 objects, red contours), from the “Gold” 
SN la sample (156 objects; blue solid lines, derived 
assuming a fixed value of Hq = 65 km s“^ Mpc“^, and 
from the subset of SNe la at 2 : > 0.9 (14 objects, blue long 
dashed lines). The three colored ellipses are the conhdence 
regions (orange: 68%; light green: 90%; yellow: 99%) for 
the combined fit of type la SN-I-GRB samples. Dashed 
lines correspond to the changing sign of the cosmic ac¬ 
celeration [i.e. q{z) = 0] at different redshifts, as labeled. 
Crosses are the centers of the corresponding contours (red: 
GRBs; blue: SNe la, white: GRB-tSN la). The black dot 
marks the ACDM cosmology. The dotted line corresponds 
to the statefinder r = 1, in this case it coincides with the 
flatness condition. 


afterglow phase, when the time scales of the prompt 
emission are irrelevant, the jet is seen as a structured 
outflow with a well defined profile. This, with some 
additional constraints on the luminosity and duration 
of the energy injection, can create a universal jet pro¬ 
file [ 13 . 

Observationally, the properties of GRB jets are not 
independent. In particular, the Ghirlanda correla¬ 
tion M tells us that if the jet opening angle (or view¬ 
ing angle) and the total energy output are known, the 
typical frequency of photons can be predicted. The 
correlation is so tight that can be used, analogous ly t o 
the stretch-luminosity correlation in type la SN |8^ . 
to standardize the GRB candle. GRBs can therefore 
be used as high redshift standard candles to constrain 
any possible evolution in the properties of the dark 
energy. 
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